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Do stochastic inhomogeneities affect dark-energy precision measurements? 
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The effect of a stochastic background of cosmological perturbations on the luminosity-redshift 
relation is computed to second order through a recently proposed covariant and gauge-invariant 
light-cone averaging procedure. The resulting expressions are free from both ultraviolet and in- 
frared divergences, implying that such perturbations cannot mimic a sizable fraction of dark energy. 
Different averages are estimated and depend on the particular function of the luminosity distance 
being averaged. The energy flux, being minimally affected by perturbations at large z, is proposed 
as the best choice for precision estimates of dark-energy parameters. Nonetheless, its irreducible 
(stochastic) variance induces statistical errors on Qa{z) typically lying in the few-percent range. 
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Establishing the existence of dark energy and deter- 
mining its parameters is one of the central issues in mod- 
ern cosmology. Evidence for a sizable dark-energy com- 
ponent in the cosmic fluid comes from different sources: 
CMB anisotropics, models of large-scale-structure forma- 
tion and, most directly, the luminosity redshift relation 
of Type la supernovae, used as standard candles. 

In this latter case, on which we concentrate our atten- 
tion, the analysis is usually made in the simplified context 
of a homogeneous and isotropic (FLRW) cosmology. The 
issue has then been raised about whether inhomogeneities 
may affect the conclusion of such a naive analysis. Inho- 
mogeneous models in which we occupy a privileged posi- 
tion in the Universe, for instance, can mimic dark energy 
(as first pointed out in pQ), but look both unrealistic and 
highly fine-tuned. More interestingly, we should address 
this question in the presence of stochastically isotropic 
and homogeneous perturbations of the kind predicted by 
inflation. We present here the main ideas and results 
of such a study, leaving to a forthcoming paper [5] their 
detailed derivation and discussion. 

There is by now general agreement that super-horizon 
perturbations cannot mimic dark-energy effects [3J. By 
contrast, the impact of sub-horizon perturbations is still 
unsettled [IH5] owing to the appearance of ultraviolet di- 
vergences while computing their "backreaction" on cer- 
tain classes of large-scale averages [6] . The possibility 
that these effects may simulate a substantial fraction of 
dark energy, or that they may at least play some role in 
the context of near- future precision cosmology, has to be 
seriously considered. 

In order to address these issues we have studied the 
luminosity-redshift relation in a spatially-flat ACDM 



model, perturbed by a stochastic background of inho- 
mogeneities. The luminosity distance d^ now depends 
on the redshift z as well as on the angular coordinates 
of the sources, and must be inserted in an appropriate 
light-cone and ensemble average |S]. Unlike for the 
cases treated in [5] E] , we find a result free of ultraviolet 
divergence and with no significant infrared contributions 
either. As a consequence, corrections are typically small, 
certainly too small to mimic a sizable fraction of dark 
energy. However, interestingly enough, both their size 
and their z-dependence strongly depend on the particu- 
lar function of d^ being averaged. 

We find, in particular, that the energy flux <f> ~ d^ 2 
is practically unaffected by inhomogeneities, while the 
most commonly used variables (like the distance modulus 
/i ~ 5 log 10 dr,) may receive much larger corrections. This 
creates (at least in principle) intrinsic ambiguities in the 
measure of the dark-energy parameters, unless the back- 
reaction of stochastic inhomogeneities is properly taken 
into account. 

Let us start by recalling the standard expression for 
the luminosity distance in an unperturbed flat ACDM 
model, with present fractions of critical density and 
O m = 1 f^A : 
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Consider now the expression for in the correspond- 
ing perturbed geometry. Combining light-cone and en- 
semble averages (denoted, respectively, by brackets and 
over-bars) , we can write the averaged result in the form: 
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where fd(z) represents the "backreaction" on due to 
inhomogeneities. For consistency, g?l has to be com- 
puted (at least) up to the second perturbative order since 
ensemble averages of first-order quantities are vanishing 
for stochastic perturbations. A detailed computation of 
fd{z) would thus enable to extract the "true" value of the 
dark-energy parameters from the measurement of {di)(z) 
after taking the correction into account. 

However, as already stressed in [8], given the covari- 
ant (light-cone) average of a perturbed (inhomogeneous) 
observable S the average of a generic function of this 
observable differs, in general, from the function of its av- 
erage, i.e. (F(S)} F((S}). Expanding the observable 
to second order as S = So + Si + S2 + • • • , one finds: 



the basis of the particular simplicity of this average and 
to its minimal deviation from the homogeneous value. 

To compute this quantity in the perturbed geometry of 
our interest, we need to express it in a gauge where the 
stochastic background of cosmological perturbations is 
explicitly known up to second order. To this purpose, we 
can use the standard Poisson gauge where we include first 
and second-order scalar perturbations, neglecting their 
tensor and vector counterparts (see [5] for a discussion of 
this point). Performing the relevant transformations to 
second order we arrive at the following analogue of ([2|: 



(d -, H(d ™ ) -, (/$(z)) _ lEE(d ™ r , [1+Mz)h 



(F(S)) = F(S ) + F'(S )(S - S ) + F"(S )(Sf/2). (3) where U has in general the following structure: 



(7) 



Thus different functions of the luminosity distance are 
differently affected by the inhomogeneities, and require 
different "subtraction" procedures. Finding the func- 
tion that minimizes the backreaction will help of course 
for a precision estimate of the cosmological parameters. 
One of the main claims of this paper is the identifica- 
tion of such an optimal observable with the energy flux 
$ = L/(And 2 L ) received from a standard candle of lu- 
minosity L located on the observer's past light-cone. We 
now illustrate how we have performed such a calculation. 

The average value of $, obviously controlled by the av- 
erage of dj 2 , has to be carried out on the past light-cone 
of the observer, at a fixed redshift z, using the gauge- 
invariant prescription introduced in [TJ. This is most 
conveniently done |5] in the so-called geodesic light- 
cone gauge, where the metric depends on six arbitrary 
functions (T,U a ,^ a b, a,b — 1,2), and the line-element 
takes the form (with 9 1 = 9, 9 2 = <•&): 

ds 2 = T 2 dw 2 -2TdwdT+ lab {de a -U a dw){d6 h ~U b dw). (4) 

In this gauge the past light-cone is defined by the condi- 
tion w = wo = const, and the redshift is given by: 



l + z = T(w Q ,To 1 a )/T(w Q ,T 1 9 a ). 



(5) 



Furthermore, the luminosity distance of the source is sim- 
ply expressed as [5] d^ = (l + z) 2 7 1 / 4 (sin6')~ 1 / 2 , yielding 
the following exact result: 
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where 7 — det 7 a f,, and t(z, 9 a ) is obtained by solving Eq. 

The above expression has a simple physical interpre- 
tation: the averaged flux, for a given z, is inversely pro- 
portional to the proper area (computed with respect to 
the metric Q ) of the surface lying on our past light-cone 
at the given value of z. Flux conservation is probably at 
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Here Zi, Zi,i, 1% are, respectively, the first-order, 
quadratic first-order, and genuine second-order contribu- 
tions of our stochastic fluctuations. After solving the rel- 
evant perturbation equations [9] they can all be expressed 
in terms of the first-order Bardeen potential ^(x, n). Us- 
ing the stochastic properties of this perturbation, and 
expanding in Fourier modes ^'^(ry), we can then obtain 
an expression for (i$) -1 where first-order contributions 
drop out because of the ensemble average, and the scalar 
perturbations only appear through the so-called dimen- 
sionless power spectrum, V{k 1 rj) = (fc 3 /27r 2 ) | (77) | 2 . 

Unfortunately, (I*) -1 contains integrals over null 
geodesies lying on the past light-cone of the given 
observer [S], which get intertwined with the time- 
dependence of V 1 forcing us to proceed with an approx- 
imate numerical integration. This will be done below, 
after inserting (as an instructive example) an illustra- 
tion of the limiting CDM case, where all integrals but 
the one over k can be done analytically thanks to the 
time-independence of V [5] . 

In that case the result can be written in the form 



/•(*) 



dk 



V(k) f 1 , 1 (k,z) + f 2 (k,z) 



(9) 



where f± t i and fi are complicated -but known- analytic 
functions of their arguments Furthermore, the lead- 
ing contribution in the region of z relevant for dark- 
energy phenomenology comes from terms of the type 
f(k,z) ~ (k/Ho) 2 f(z), where Ho is the present Hubble 
scale. We can then write, to a very good accuracy, 



U(z)~ fx,x{z) + h{z) 
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FIG. 1: The fractional correction /$ of Eq. (|9| (solid curve), 
compared with the same quantity given to leading order by 
Eq. (JXoJ) (dashed curve), in the context of an inhomogeneous 
CDM model. We have used for V(k) the inflationary scalar 
spectrum with the WMAP parameters [10] and the transfer 
function given in [11] (see also [8]). The plotted curve refers, 
as an illustrative example, to an UV cutoff kuv = IMpc -1 . 



FIG. 2: The fractional correction to the flux /$ of Eq. Q 
(thin curves) is compared with the fractional correction to the 
luminosity distance fd of Eq. ( |13[ ) (thick curves), for a ACDM 
model with Q\ = 0.73. We have used two different cutoff 
values: kuv = O.lMpc -1 (dashed curves) and kuv = IMpc -1 
(solid curves); the spectrum is the same as that of Fig. [TJ 
adapted to ACDM. 



where an explicit calculation gives [2J: 
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The absolute value (and sign) of /$(z) are illustrated in 
Fig. [TJ showing the accuracy of the leading order terms 
( [To] ) and confirming that the backreaction of a realistic 
spectrum of stochastic perturbations induces negligible 
corrections to the averaged flux at large z (the larger 
corrections at small z, due to "Doppler terms", has been 
discussed in [8]). In addition, it shows that, in any case, 
such corrections have the wrong z-dependence (in par- 
ticular change sign at some z) to simulate even a tiny 
dark-energy component. For the considered spectrum 
(behaving as k ns ~ 5 log 2 k at large k, see [TT]) the spec- 
tral integral is convergent and very weakly sensitive to 
the chosen value of the UV cutoff [8] representing here 
the limit of validity of our perturbative approach. 

We now come to the more realistic ACDM case, where 
the /$ correction should be obtained by a full numer- 
ical integration of Eqs. (7[8[. For simplicity, we will 
only take into account those terms giving the leading (k 2 - 
enhanced) contributions in the CDM case. For ACDM 
we can generally expect a smaller correction due to the 
fact that the spectrum is now suppressed, at large k, by 
a lower value of the equality scale k eq [TT] . This is con- 
firmed by the explicit numerical result for |/$| presented 
in Fig. [2j The small value of |/$| at large z leads us 
to conclude that the averaged flux is a particularly ap- 
propriate quantity for extracting from the observational 
data the "true" cosmological parameters. As we are go- 



ing to see now, the situation is somewhat different for 
other functions of d^. 

Indeed, let's apply the general result ^ to the flux 
variable, S = $, and consider two important examples: 
F($) = $ -1 / 2 - d L , and = -2.5 log 10 $ + const - 

fi (the distance modulus). For the luminosity distance, 
following the notations of Eq. ^ and using the general 
result ([3]), we obtain: 



/ d = -(l/2)/ $ + (3/8)((<f> 1 /$ ) 2 ). 
Similarly, for the distance modulus we obtain: 
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= -1.25(log 10 e) 2/» -(($1/$ 
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(14) 



where is defined in Eq. Q. 

As clearly shown by the two above equations, the cor- 
rections to the averaged values of rfi and \x are quali- 
tatively different from those of the flux (represented by 
/$), because of the extra contribution (inevitable for any 
non- linear function of the flux) proportional to the square 
of the first-order fluctuations. As mentioned before, the 
averaged flux corrections have leading spectral contribu- 
tions of the type k 2 V(k)\ on the contrary, the new cor- 
rections to 4 and fi are due to the so-called "lensing 
effect", dominate at large z, and have leading spectral 
contributions of the type k 3 V(k) (as already discussed 
in [5]). The explicit numerical integration, reported in 
Fig. 2, confirms that, as a result, |/$| -C fd at large z. 
We stress that even the /c 3 -enhanced contributions are 
UV-finite for the case under consideration. 

Let us now briefly discuss to what extent the enhanced 
corrections due to the squared first-order fluctuations 
can affect the determination of the dark-energy param- 
eters if quantities other than the flux are used in the 



4 



fits. To this purpose we consider the much used (av- 
erage of the) distance modulus given in Eq. (14), re- 



ferred as usual to the homogeneous Milne model with 
fi M = 51og 10 [(2 + z)z/(2H Q )]. In Fig. [|we compare the 
averaged value (/i) — fi M with the corresponding expres- 
sion in a homogeneous ACDM model with different val- 
ues of J7a- We also show the expected dispersion around 
the averaged result, represented by the square root of the 
variance [8]. The latter is given by: 



<M 2 )-(MJ =±2.5(log 10 e)V(($i/$o)0; (15) 
while for the flux we simply find: 

y (($/$ ) 2 ) - (w^>) 2 = iv^i/^o) 2 ) . (i6) 

As illustrated in Fig. 3, we find that, even for the 
distance modulus, the effect of inhomogeneities on the 
average only affects the determination of £Ia at the third 
decimal figure (see also Fig. 2 and [2] ) , at least for the in- 
flationary power spectrum with the ACDM transfer func- 
tion of [TT]: in that case, the curves for (fj) and ^t FLRW 
are practically coincident at large z. We have considered 
other spectra which take into account non-linear effects 
and have more power at short scales, like those obtained 
following [12] . Using such spectra only affects very mildly 
the fc 2 -enhanced terms (hence the flux) while they in- 
crease the corrections wherever the fc 3 -enhanced lensing 
terms play a major role. In particular, the variance due 
to the fluctuations, which is already at the few-% level at 
large z for the power spectrum of [H] (see Fig. [3]) , can be 
further increased [2]. Note that, even for these improved 
spectra, all our integrals are still free of UV divergences 
since, in any case, V falls faster that fc -3 (i.e. the matter 
density constrast spectrum grows slower than k). 

Our main conclusions can be summarized as follows: 

(1) Unlike in the case of other large-scale averages 
|B] , a physically motivated light-cone averaging of the 
luminosity-redshift relation leads to expressions for the 
backreaction of stochastic inhomogeneities that are free 
from UV (as well as IR) divergences. This allows us to 
conclude that such perturbations cannot simulate a sub- 
stantial fraction of dark energy. 

(2) The actual value of the backreaction strongly de- 
pends on the quantity being averaged. It turns out to be 
minimal for the flux <E>, which is also practically insensi- 
tive to the short-distance behaviour of the power spec- 
trum. Therefore, the flux stands out as the safest observ- 
able for precision cosmology. For other observables, such 
as the distance modulus, the backreaction is considerably 
larger and is more sensitive to the spectrum used. 

(3) The dispersion due to stochastic fluctuations is 
much larger than the backreaction itself, implying an ir- 
reducible scatter of the data that may limit to the percent 
level (see Fig. [3]) the precision attainable on cosmological 
parameters because of the present limited statistics. 
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FIG. 3: The averaged distance modulus {/i} — fx (thick solid 



curve), and its dispersion of Eq. (15 1 (shaded region) are 
computed for Qa = 0.73 and compared with the homogeneous 
value for the unperturbed ACDM models with Qa = 0.69, 
0.71, 0.73, 0.75, 0.77 (dashed curves). We have used k uv = 
IMpc -1 and the same spectrum as in Fig. [2] 



(4) We calculated here the full second order effect of 
stochastic perturbations. Possible contributions coming 
from the non-perturbative regime on length scales much 
smaller than 1 Mpc have still to be taken in consideration 
before final conclusions can be drawn. 
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